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O ' Abstract 






We study the infrared part of the spectrum for UV-filtered staggered Dirac operators and 

QJ^- compare them to the overlap counterpart. With sufficient filtering and at small enough 

^^ . lattice spacing the staggered spectra manage to "mimic" the overlap version. They show a 

^^ I 4-fold near-degeneracy, and a clear separation between would-be zero modes and non-zero 

modes. This suggests an approximate index theorem for filtered staggered fermions and 

^ ■ a correct sensitivity to the topology of QCD. Moreover, it supports square-rooting the 

K* i staggered determinant to obtain dynamical ensembles with Nf = 2. 

(N 

^ ; 1 Introduction 

o 

Q . Recent years have brought a remarkable increase of interest in dynamical {Nf = 2 or Nf = 2+1] 
QCD runs employing staggered sea-quarks. One of the conceptual issues with these efforts 
is that the staggered quark action is not doubler-free; a single staggered field will eventually 



c^ 






p '_ generate four flavors in the continuum. Hence, to create an ensemble with Nf = 2, one simply 
uses the square-root of the staggered determinant, and for A^ = 2-1-1 an additional quartic-root 
with a higher mass. The problem is that taking a fractional power of the determinant is, in 
general, not a legitimate operation in quantum fleld-theory, contrary to positive integer powers. 
rN I The worst-case scenario is, therefore, that these runs might not represent an ab-initio approach 

o3 I to QCD in the strong coupling regime. In this respect it does not help that the recent staggered 
A^ = 2-1-1 results with light sea-quarks look very promising jT] . It has been clearly pointed out by 
Jansen at the Tsukuba conference j^ and in DeGrand's review jH] that - if the square-rooting 
trick cannot be justifled - dynamical staggered runs with Nf^AlN must be considered a model 
of QCD. The results in P show that it would be a much better model than quenched QCD, 
but from a fundamental point of view little is won if the answer would be that there is no 
justification for the fractional power. 

The goal of this note is to confront the low-lying eigenvalues of the massless staggered 
Dirac operator j4j with analogous spectra of a "benchmark" operator which is known to be 
without conceptual problems. We choose the massless overlap operator [Hj, since it satisfies 
(like domain-wall ^ and classically perfect [Zj fermions) the Ginsparg- Wilson relation jHj 

D^, + ^,D = -D^,D (1) 

P 



with p a parameter to be specified later. This relation is of prime interest, because it implies 
invariance of the action (at finite lattice spacing) under the adapted chiral transformation Pj 

5ilj = -i,{l--D)ij, # = ?^75, (2) 

P 

which in turn excludes additive mass renormalization and prevents operators in different chiral 

multiplets from mixing. On the other hand, the staggered action breaks the full SU{Nf = 4)A 

flavor symmetry, but the remnant ?7(l)-symmetry 

still protects the fermion mass against additive renormalization. 

A point worth emphasizing is that modern staggered simulations use highly "improved" 
staggered quarks. These are operators where one replaces in the usual definition 

^stag = Wr^,ix){u,ix)6,^f,,y - Ulix-f,)5,.f,,y) (4) 

[with 77^ (x) = (— 1)^'^<''^'^] the parallel transporter U^{x) by a weighted sum of several gauge- 
covariant paths from x to x + jl, and there is an abundance of proposals which terms should 
enter and what are useful weights JU]- From the Symanzik point of view [llj, these actions are 
in general not improved - they are in the same class with 0{a^) cut-off effects as the original 
staggered action (albeit the extrapolation slope is typically reduced, and the scaling window 
might begin earlier). Since these modifications are designed to reduce the dramatic background 
fluctuations at the cut-off level (in particular the "taste-changing" interactions due to highly 
virtual gluon exchanges [IHl), one should rather speak of "UV-filtered" staggered quarks. 
The massless overlap operator is defined through |SJ 

D— = p(l + D^^pwt^w )-i/2) ^ ^(1 ^ ^^ sign(75l^^,)) (5) 

with D"^ the Wilson operator at negative mass —p. In the free field limit p=l is appropriate, 
but for finite (3 one usually chooses some p > 1 to optimize the locality of _D°™'' p2j. Given 
the success of the UV-filtered staggered operators, it is natural to consider modified overlap 
operators where the Wilson Kernel is constructed from smoothed parallel-transporters, too [l3l 
ITi] . The resulting filtered overlap operator is still in the Symanzik class with Ola"^) artefacts, 
if one stays at a fixed smearing level for all /?. In other words: The choice of the covariant 
derivative in D"^ is an O^a^) ambiguity in _D°™'' like the choice of p. We thus take the liberty 
to fix p = l at all /5- values (except for one check), and to just change the smearing level. 

With this setup we are ready for an investigation, which mainly follows the "spectral hint" 
section of Ref. ^3]. To come back to the starting point, what one would like to check is whether 
UV-filtered staggered quarks develop approximately 4-fold degenerate eigenvalues which coin- 
cide with the overlap counterparts such that the rooted staggered determinant generates, up 
to Oi^a"^) cut-off effects, the correct ensemble with 1 or 2 light sea-quarks, 

det D"™"^ oc (det D^'^s) V4 + Qia^) . (6) 

Our approach is to focus on the technically simpler version 

det I^&T-mtcred oc (det D'^^^^.^^J'/^ + 0{a') , (7) 

which is equivalent, since the two left- and the right-hand sides differ by a trivial rescaling and 
O(a^) terms only. If such a relationship can be shown, the overlap operator would be a local 
operator with the same determinant, up to cut-off effects, as the rooted staggered operator, 
hence providing a bypass to the locality issue for the valence quarks presented in 



2 Low-lying eigenvalues of UV-filtered Dirac operators 

The plan is to compare - configuration by configuration - the low-lying eigenvalue spectrum of 
a UV-filtered staggered Dirac operator against that of a filtered overlap operator. 

Regarding the filtering, we decided not to follow the details of current state-of-the-art stag- 
gered fermions, but to concentrate on two prototype constructs, the APE [in] and HYP |T7j 
smoothed hnks with the S'[/(3)-projection of Ref. [IH] which are used for both types of lattice 
fermions. Throughout, we used the parameters q;ape = 0.5 [12] and aHYP = (0.75, 0.6, 0.3) [T7j . 
By comparing the two recipes and the various iteration levels, we will be able to judge how 
much our results are robust against changing the details of the smearing procedure. 

The key issue is whether these filtered staggered Dirac operators develop an approximate 
4-fold degeneracy, as it has been demonstrated in great detail in 2D (where the degeneracy is 
2-fold) Pl], and seen more recently in 4D [THj . 

2.1 Dynamical lattices 

Let us begin with a full QCD set that was used for spectroscopy. CU_0001 is an ensemble of 
49 configurations with 16'^ x 32 geometry, generated (via the square-root trick) with 2 "tastes" 
of staggered quarks at bare mass m = 0.01 and coupling (3 = 5.7. The lattice spacing is around 
a~0.1fm. We downloaded it from [^U], the pertinent original publication is ^TD. 

In order to select configurations with and without overlap zero modes, we first measured a 
lattice version of the continuum topological charge gnai = -^^J2x'^''^{^fj,uapFfj,uF^p), where F^,^ 
is defined via the clover-leaf construction involving smeared links (7 HYP steps). Based on 
this, we selected the configurations shown in Tab. ^ for the spectral analysis. The occurrence 
of configurations where the overlap index n_—n^ = ^tr(75D°^'^'^) (denoted q in the following) 
depends on the smearing level and type (or on the projection parameter p) is of no concern - it 
is a typical 0{a^) effect. The ambiguity is not removed by using vigorously filtered operators, 
but it disappears in the continuum limit. 

Since D^^'^^ and Z}°^'='' are normal, we determined the spectra of the hermitean operators 
jjsta.gi j~)sta.g ^ jjover't jjover ^^^^ recoustructcd the Staggered and overlap eigenvalues, respectively. 
The eigenvalues of the hermitean operators have been determined with the Ritz functional 
method [22] . We present the eigenvalues of D^^^^ and the chirally rotated 



A=(l/A-l/(2p))"' 



of 1}°™^" which are both purely imaginary. Figs. ^121121 contain our results. The thin-link 
operators are compared on the left-most side, and the filtered cousins with the same type 
of parallel-transporters are arranged next to each other. With a~^ ~ 2 GeV the energy axis 
includes modes up to 160 MeV. 
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Table 1: Original name, nickname, plaquette, gnai and n_ — n+ of the CU-0001 configurations. 
Whenever they would not agree the indices for the standard, 1 APE, 3 APE, 1 HYP and 3 HYP 
overlap operators are given separately. 
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Figure 1: Eigenvalues on the configurations 00 and 06 in the ensemble CU_0001. 
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Figure 2: Eigenvalues on the configurations 14 and 26 in the ensemble CU_0001. 
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Figure 3: Eigenvalues on the configuration 38 and 48 in the ensemble CU_0001. 
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The original operators show no similarity at all. Upon employing one APE-smearing step, 
the staggered spectra get stretched while the overlap-spectra get squeezed. Still, there is no 
similarity, but the eigenvalue densities (divided by 4 in the staggered case) have been equalized 
quite a bit. At intermediate filtering levels (3 APE and 1 HYP steps) all staggered eigenvalues 
have grouped into almost-degenerate twin-pairs. At the 3 HYP level, finally, a second grouping 
of twin-pairs into quadruples of near- degenerate staggered eigenmodes has taken place. Stretch- 
ing the staggered 3 HYP spectra with a factor Z~1.2, one finds a rather convincing agreement 
between the quadruples and the (individual) overlap modes. There is a clear gap between 
staggered "would-be" zero modes and non-zero modes, though this rule is challenged by some 
configurations. On conf_014 (and to some extent on conf_038) the operator -Dshyp reproduces 
the confusion that the overlap has regarding the charge of this configuration (cf. footnote ^). 
And conf_048 seems to be hard to resolve, simply for its high |g|. Note, finally, the stability 
(on a clear majority of configurations) of the low-lying overlap spectrum from 3 APE through 
3 HYP. 

All these findings are in complete analogy to what has been seen in great detail in 2D [T5I . 
The only difference is that the near-degeneracy is 2-fold there, since Q is the lattice remnant 
of the SU{Nf = 2) symmetry in 2D. 

2.2 Matched quenched lattices 

To check whether O(a^) effects really capture the difference between filtered staggered and 
overlap spectra, we generated 4 sets of matched quenched lattices, aiming at L/ro = 2.24 with 
the Sommer scale tq ~ 0.5 fm |2S1- In a physical volume V ~ (1.12 fm)^ ~ 1.57fm there are 
finite volume effects, but we expect them not to spoil the spectral analogy, if it emerges at the 
given coupling. Using the interpolation formula in |21] for the Wilson gauge action, one finds 
that the pairs (L = 6,/? = 5.66), (L = 8,/5 = 5.79), (L = 12,/3 = 6.0), (L = 16,/? = 6. 18) produce the 
desired volume (the first coupling is slightly out of bound, the formula holds for 5.7 </3< 6.57). 

Tables |2l-El contain the absolute values of the overlap indices on our 6^, 8^, 12^, 16^ lattices. 
As expected, the fraction of configurations on which the 1 APE, 3 APE, 1 HYP and 3 HYP 
overlap operators agree on the charge increases with (3 (for j3 — > cxd, this fraction is one). 
TableElis interesting in another respect: The normal "rule" (at high j3) that additional smearing 
operations will reduce the charge (here, we ignore the unsmeared column) is broken at low (3 
(at 13 = 5.66 there are five such configurations). Finally, Table El also contains our test regarding 
the impact of the projection parameter p in (jHl). On the first configuration, besides using p=l, 
we determined the spectrum with p = 1.6. The standard and 1 APE overlap operators prove 
sensitive to this shift, the more vigorously filtered versions not. This is fully compatible with 
the view that the charge ambiguity is an Ola"^) effect; with strong UV-filtering the coefficient 
in front of the artefact gets smaller. Interestingly, -DgsAPE 63 hyp ^^i^^ yield reasonable charge 
determinations, in spite of these operators being entirely non-local on our small lattices. 

We tested on the 8^ configurations whether the S'f/(3)-projection in the APE ^Hl or 
HYP ^7j smearing recipes could be skipped. The answer is negative: With one smearing 
step, the impact is not so dramatic, but after two more the staggered spectra are dramati- 
cally squeezed (rather than stretched out), so that there is no agreement with the (filtered and 
projected) overlap spectrum at all. 

Contrary to the CU_0001 case, we refrain from showing all spectra. For the lack of a good 
criterion, we simply select the first configuration of our 6^, 8^, 12^, 16^ chains, respectively. 
In Figures El and jSl the (ordered) staggered eigenvalues (with projection) are plotted against 
their serial number. One clearly observes how the expected 4-fold near-degeneracy gets more 
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Table 2: Number of zero modes of various overlap operators on our quenched 6^ configurations. 
The 1 APE, 3 APE, 1 HYP and 3 HYP operators agree for 4 of the 10 conhgurations. 
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Table 3: Number of zero modes of various overlap operators on our quenched 8^ configurations. 
On the hrst configuration, the values for p = 1.6 are in brackets. The lAPE, 3 APE, IHYP 
and 3 HYP operators agree for 5 of the 6 configurations. 
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Table 4: Number of zero modes of various overlap operators on our quenched 12"^ configurations. 
The 1 APE, 3 APE, 1 HYP and 3 HYP operators agree for all 6 configurations. 
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Table 5: Number of zero modes of various overlap operators on our quenched 16^ configurations. 
The 1 APE, 3 APE, IHYP and 3 HYP operators agree for all 4 configurations on which they 
have all been determined. 
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Figure 4: Staggered eigenvalues with/without relative normalization factor on the first one of 

our 6'^ and 8^ lattices, compared to the corresponding overlap spectrum. 
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pronounced as (3 increases; at /3 = 6.18 there is a clear "terrace dynamics" . This picture holds for 
-^iHYP) ^^d fo'^ -^3 HYP i^ is even more convincing. On the left, the IR-spectrum of the overlap 
operator with the same type of filtering is shown. Finally, the appropriate relative normalization 
factor (slope in Figs. IHl and [71 see below) is applied to the staggered spectrum to make it a fair 
comparison. The bottom line is that, at sufficiently weak coupling, it is appropriate to identify 
the geometric mean of a near- degenerate quadruple with a single eigenvalue of the overlap 
operator. Whether there is, in addition, a relationship between the modes - ideally, one would 
like to give a construction which, in analogy to the one by Kogut and Susskind, manages to 
(approximately) thin out the degrees of freedom - is a challenge for the future. 

Our Figures |3] and El cover a fixed energy range up to ~ 370 MeV. They demonstrate that 
both the separation between would-be zero modes and non-zero modes and the grouping into 
quadruples crucially depend on the lattice spacing - below a^^ ~ 2 GeV even a high level of 
filtering does not help. And, by comparing to the filtered overlap we see a quantitative agreement 
for a-^>2-2.5GeV. 

2.3 Correlation between staggered and overlap eigenvalues 

To check whether the quartic-rooted filtered staggered operator really reproduces the spectrum 
of the filtered overlap one may plot them, eigenvalue by eigenvalue, against each other. 

Restricting ourselves to the positive half of the spectrum, we group the staggered eigenvalues 
into two classes: The would-be zero modes (we name them 0<Ci<C2<---) and the (true) non- 
zero modes (denoted by Ai < A2 < ... with Cmax=2|g| < Ai). For the n-th would-be zero mode 
we take A/C2n-iC2n, since (— C2, — Ci) Ci? C2) is the substitute of the first zero mode. The m-th 
non-zero mode is associated with (A4m-3 ■ ••• ■ ^imY^^- This classification is a strict function of 
q determined via the overlap against which one compares. 

These geometric means are then plotted against the overlap eigenvalues A, results being 
shown in Figs. EHZl Here, every symbol corresponds to an individual configuration. 

What one sees immediately, is that the agreement gets much better with increasing (3. In 
three graphs there is a "pile" of would-be zero modes above the origin. These are classified so, 
because the overlap finds a non-zero index, but then they happen to be not-so-small. Since the 
overlap operator itself has an 0{a^) ambiguity^, this phenomenon is expected to become rare 
at weak coupling. Unfortunately, even at /5 = 6.18 we have one configuration with such a charge 
ambiguity. The pile also shows up when one overlap is plotted versus another one, say -Dihyp 
vs. -DyHYP- O^ ^^6 other hand, the fact that one configuration in the j3 = 6.0 ensemble goes 
astray is due to UV-noise which does not affect the charge. One learns that at this coupling 
even the 1 HYP operator suffers from large UV-fluctuations in its IR-spectrum. 

Switching in Figs. IHl El from the 1 HYP level (black) to 7 HYP staggered vs. 7 HYP overlap 
(light), the correlation gets tighter and the relative normalization factor Z moves closer to 1. 
Increasing, on the staggered side, the smearing level to 7 HYP, but staying, on the overlap 
side, at IHYP (i.e. D^-^yp ^^- -^ihyp)' the Z-factor still moves towards 1, but the correlation 
does not improve (not shown). We plan to check that the tight correlation between the lowest 
eigenmodes extends to the whole determinant, i.e. to verify eqns. (0121) directly. 

Once f5 is so large that the pile of not-so-small would-be zero modes is gone and a clear 
gap to the non-zero modes is observed (cf. Fig. EJ, the filtered staggered fermions will satisfy 
an approximate index theorem. This has been demonstrated in great detail in 2D [14 , and 
Refs. (23 Uni have added evidence that this holds in 4D, too. 

^With another p or an alternative smearing level a different number of staggered modes might have been 
identified as would-be zero modes. 
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Figure 6: Quartic-rooted staggered versus overlap modes at /3 = 5.66 and j3 = 5.79. Different 

symbols refer to different configurations. 
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Last but not least, the mode- by-mode correspondence implies that (filtered) staggered 
fermions on sufficiently fine lattices reach agreement with the prediction of random matrix 
theory, as it has been observed for overlap fermions J2E]- 

2.4 Scaling analysis for eigenvalue-splitting pseudo-observables 

To define the fuzziness of the lowest would-be zero mode (if present) we form 






To define the fuzziness of the first (true) non-zero mode we construct 



A _ \/A4A3 — VA2A1 , 

' ~ (A8A7A6A5)V4-(A4A3A2Ai)V4- ^ ^ 

These expressions quantify the non-degeneracy of a quadruple on the scale of its distance to 
the next one, and tend to zero if the four-fold degeneracy becomes exact. 

If the non-degeneracy of the staggered eigenmodes and the fuzziness of the would-be zero 
modes are just 0{a?) artefacts (though, numerically, they may be so dramatic that without 
UV-filtering not even the faintest analogy is seen at accessible couplings), one would expect 
that the pseudo-observables ()10|) and Q vanish (asymptotically) in proportion to a?. 

Fig. IHl contains our results for the 8^, 12^, and 16"^ lattices (the points for the 6^ geometry 
were omitted since they are definitely not in the scaling regime). Given the smallness of our 
samples, we have opted for the median instead of the arithmetic mean - this estimator is 
biased, but far more robust against outliers. Indeed, the 1 HYP and 3 HYP data seem to be 
consistent with the hypothesis of asymptotic scaling - although our data do, of course, not 
pin down the exponent. A linear fit (versus a^) through the L = 16, 12 points (constrained to 
go through 0, i.e. with 1 d.o.f.) yields an acceptable x^y and its extension to the right seems 
compatible with the L = 8 data. On the other hand, the unfiltered data are essentially fiat, 
and there is no way to fit them in this manner. This is precisely the scenario mentioned in the 
introduction. While we have no doubts that - moving to much weaker couplings - even the 
unsmeared staggered operator's pseudo-observables Ao,i would eventually turn down and reach 
the origin linearly, the big difference with the 1 HYP and 3 HYP versions is that this happens 
at accessible couplings. 



Furthermore, the lesson from 2D is that a typical JC,2qC,2q-i of the massless operator gives 
a reliable estimate down to which quark mass the massive staggered Dirac operator yields 
trustworthy results [14j. 

2.5 Technical bonus 

It is interesting to note that the construction of the UV-filtered overlap operator is much 
cheaper than the standard version. For the unsmeared CU configurations the condition number 
of H^ = (75/^^^)^ was in the range 10^ — 10®. After projecting out the 11 lowest eigenmodes, 
this figure gets reduced to ~350 — 450, and the degree of the associate Chebychev polynomial 
is ~ 400. After 3 HYP steps, the condition number is ~140 — IS'OOO and ~ 60 — 100, before 
and after projection, respectively, and the required degree of the Chebychev polynomial on the 
subspace orthogonal to the 11 lowest modes is ~60. 

Note that this impact on the degree of the polynomial is seen after the lowest 11 eigenvectors 
have been projected out. This strengthens the impression that one gets from comparing Fig. 2 
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in Ref. [22] to Fig. 3 in Ref. [2H]. In the latter case, the standard (unfiltered) Wilson kernel is 
evaluated on pure gauge backgrounds generated with the Iwasaki/DBW2 action, and the lowest 
eigenvalues of H^ get considerably lifted compared to the Wilson gauge action case. However, 
after a fixed number of eigenmodes have been projected out, this advantage is considerably 
reduced (see Fig. 3 of Ref. |2HI), since the bulk part of the spectrum is less sensitive to the 
improvement coefficient than the bottom end of the eigenvalue distribution. In Ref. J22I the 
fermion action is changed, and this has a huge impact even after a fixed number of eigenvectors 
is projected out. In line with this, we observe that the hulk of the eigenmodes of H^ for a HYP 
smeared gauge configuration (which we see as a filtered fermion on the original background) 
starts at significantly higher values. This holds, regardless whether this configuration was 
produced with the Wilson gauge action or stems from a full QCD ensemble. Therefore, the 
construction of the overlap operator with a UV-filtered kernel is considerably cheaper than for 
a standard Wilson kernel even after projecting out the few lowest eigenmodes of H^. 

3 Summary and Outlook 

We have studied a class of UV-filtered staggered and overlap operators that are built with 
standard techniques, comparing the IR-part of their spectra with and without filtering. We 
close with a summary of our findings and a few comments on their implications. 

1. With sufficient filtering and on fine enough lattices, the staggered fermions develop an 
ability to separate would-be zero modes from non-zero modes, and arrange them in near- 
degenerate groups of 4. The second point is that both the number of zero-modes and the 
quadruples reproduce, up to 0{a?) effects, the findings of the filtered overlap operator. 
This implies an approximate index theorem for filtered staggered quarks. 

2. For both aspects, the "critical" coupling is (in quenched Wilson terms) around /9w = 6.0. 
Below, the agreement vanishes rapidly, and there is no way to compensate for this by 
resorting to higher smearing levels. 

3. The agreement, at the level of a single configuration, with the low-energy spectrum of the 
overlap operator implies that UV-filtered staggered fermions will eventually reproduce 
the correct random-matrix-theory universality class and the correct continuum limit of 
topological susceptibility and associate fiavor-singlet observables, as a corollary. 

4. A rudimentary scaling analysis of two pseudo-observables designed to capture the fuzzi- 
ness of the staggered would-be zero modes and non-zero modes, respectively, suggests 
that our data are consistent with the hypothesis that these staggered artefacts vanish 
like 0{a?). This, together with the quantitative agreement with the corresponding over- 
lap eigenvalues, forms evidence in favor of using the fractional power ((Zj) of the filtered 
staggered determinant. 

5. An observation important for quenched overlap spectroscopy is that the filtered overlap 
is much cheaper than the unfiltered one. This is a distinctive feature of the kernel in (jHI), 
and not tied to any particular action used to generate the background. 

6. Our line-up and correlation plots at weak coupling (Figs. ElandEj) bear the promise to 
trigger a new development in lattice QCD. With such an excellent one-to-four agreement 
of the low-lying spectra, it might be possible to run a dynamical Nf = 2 staggered simu- 
lation, performing the final accept/reject-step with the overlap action (or, alternatively, 
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to reweight the ensemble to A^ = 2 overlap flavors). In other words, the perspective is to 
perform a fully fledged dynamical overlap or domain-wall study essentially at the costs 
of a staggered run. 

Acknowledgments 

S.D. and U.W. would like to acknowledge stimulating discussions with Karl Jansen, Ch.H. with 
Laurent Lellouch and Leonardo Giusti. S.D. is supported by DFG in SFB/TR-9, Ch.H. is 
supported by EU grant HPMF-CT-2001-01468. 

References 

[1] C.T.H. Davies et al. [HPQCD Collaboration] , |hep-lat/0304004[ 

[2] K. Jansen, Nucl. Phys. Proc. Suppl. 129, 3 (2004) Ihep-lat/OailOSQ] . 

[3] T. DeGrand, Int. J. Mod. Phys. A 19, 1337 (2004) |hep-ph/031224l] . 

[4] L. Susskind, Phys. Rev. D 16, 3031 (1977). 



[5] R. Narayanan and H. Neuberger, Nucl. Phys. B 412, 574 (1994) |hep-lat/9307006 . 
R. Narayanan and H. Neuberger, Nucl. Phys. B 443, 305 (1995) hep-th/9411108 . H. Neu- 
berger, Phys. Lett. B 417, 141 (1998) |hep-lat/9707022 . H. Neuberger, Phys. Lett. B 427, 
353 (1998) hep-lat/9801031 . 

[6] D.B. Kaplan, Phys. Lett. B 288, 342 (1992) |hep-lat/920601^i| . Y. Shamir, Nucl. Phys. B 
406, 90 (1993) |hep-lat/9303005j. 

[7] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B 414, 785 (1994) |hep-lat/9308004 . T. De- 
Grand, A. Hasenfratz, P. Hasenfratz and F. Niedermayer, Nucl. Phys. B 454, 587 (1995) 
hep-lat/9506030 . U.J. Wiese, Phys. Lett. B 315, 417 (1993) hep-lat/93 06003j . W. Bie- 
tenholz and U.J. Wiese, Nucl. Phys. B 464, 319 (1996) |hep-la t/9510026 . 

[8] P.H. Ginsparg and K.G. Wilson, Phys. Rev. D 25, 2649 (1982). 

[9] M. Llischer, Phys. Lett. B 428, 342 (1998) |hep-lat/9802Qll|. 



[10] T. Blum et a/., Phys. Rev. D 55, 1133 (1997) lhep -lat/960 9036 . G.P. Lepage, Phys. Rev. 
D 59, 074502 (1999) |hep-lat/9809157 . K. Orgi nos, D. To ussaint and R.L. Sugar [MILC 
Collaboration], Phys. Rev. D 60, 054503 (1999) |hep-lat/9 903032 . 

[11] K. Symanzik, Nucl. Phys. B 226, 187 (1983). K. Symanzik, Nucl. Phys. B 226, 205 (1983). 

[12] P. Hernandez, K. Jansen and M. Liischer, Nucl. Phys. B 552, 363 (1999) |hep-lat/9808010| . 



[13] T. DeGrand [MILC collaboration], Phys. Rev. D 63, 034503 (2001) plep-lat/0007046| 
[14] S. Diirr and Ch. Hoelbling, Phys. Rev. D 69, 034503 (2004) |hep-lat/0311002] . 
[15] B. Bunk, M. Delia Morte, K. Jansen and F. Knechth, |hep-lat/0403022l 
[16] M. Albanese et al. [APE Collaboration], Phys. Lett. B 192, 163 (1987). 



17 



[17 

[is: 

[19 
[20 

[21 
[22 
[23 
[24 

[25 

[26 
[27 



A. Hasenfratz and F. Knechtli, Pliys. Rev. D 64, 034504 (2001) fhep-lat/0103029l. 
A. Hasenfratz, R . Hoffmann and F. Knechtli, Nucl. Phys. Proc. Suppl. 106, 418 (2002) 
|hep-lat/0110l68| . 

J.E. Kiskis, R. Narayanan and H. Neuberger, Phys. Lett. B 574, 65 (2003) 
pIep4at/0308033 . 



E. Follana, A. Hart and C.T.H. Davies, jhep-lat/0406010 



The Gauge Connection at the US National Energy Research Scientific Computing Center, 
|http : //qcd . nersc . gov/ 



F.R. Brown et al, Phys. Rev. Lett. 67, 1062 (1991). 

T. Kalkreuter and H. Simma, Comput. Phys. Commun. 93, 33 (1996) |hep-lat/9507023| . 

R. Sommer, Nucl Phys. B 411, 839 (1994) |hep-lat/ 93 100221 . 

M. Gu agnelh, R. Sommer and H. Wittig [ALPHA collaboration], Nucl. Phys. B 535, 389 
(1998) |hep-lat/9 806005|. 

P.H. Damgaard, U.M. Heller, R. Niclasen and K. Rummukainen, Phys. Rev. D 61, 014501 
(2000) |hep-lat/9907019| . 



L. Giusti, M. Llischer, P. Weisz and H. Wittig, JHEP 0311, 023 (2003) fhep-lat/0309"T89 



P. Hasenfratz, S. Hauswirth, T. Jorg, F. Niedermayer and K. Holland, Nucl. Phys. B 643, 
280 (2002) |hep-lat/0205010| . 



[28] K. Jansen and K.-I. Nagai, JHEP 0312 (2003) 038 |hep-lat/030 5009l. 



18 



